Weak coupling study of decoherence of a qubit in disordered magnetic environments 
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We study the decoherence of a qubit weakly coupled to frustrated spinbaths. We focus on spin- 
baths described by the classical Ising spin glass and the quantum random transverse Ising model 
which are known to have complex thermodynamic phase diagrams as a function of an external 
magnetic field and temperature. Using a combination of numerical and analytical methods, we show 
that for baths initally in thermal equilibrium, the resulting decoherence is highly sensitive to the 
nature of the coupling to the environment and is qualitatively different in different parts of the phase 
diagram. We find an unexpected strong non-Markovian decay of the coherence when the random 
transverse Ising model bath is prepared in an initial state characterized by a finite temperature 
paramagnet. This is contrary to the usual case of exponential decay (Markovian) expected for 
spin baths in finite temperature paramagnetic phases, thereby illustrating the importance of the 
underlying non-trivial dynamics of interacting quantum spinbaths. 
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I. INTRODUCTION 

The understanding and control of the decoherence of 
small quantum systems is central to a lot of recent de- 
velopments in the fields of nanotechnology and quantum 
computersi"— . For example, the efficacy of qubits, the 
basic building block of a quantum computer, which can 
be spin qubits, Josephson junction qubits, charge or flux 
qubits depends largely on the environments, or baths, to 
which they are often weakly coupled. The effect of envi- 
ronments on the coherence of the qubit has been studied 
in various contexts with particular emphasis on bosonic 
baths^— . In the past few years, the realization of solid 
state qubits in semiconducting heterojunctions has also 
resulted in the study of spin baths constituted of spins. 
At very low temperatures, the decoherence engendered 
by a spin environment is expected to dominate the loss 
of the coherence arising from a coupling to phononic de- 
grees of freedom. 

In the limit of weak coupling between the qubit and 
the bath, the effect of intrabath interactions have been 
explored earlier in other works^"— . In some cases, inter- 
actions were found to decrease the rate of decoherence, 
though the exact opposite was seen in cases where the 
spin bath was on the verge of a standard magnetic quan- 
tum phase transitiorti^. On the other hand, the qualita- 
tive nature of the decoherence i.e., whether it is Marko- 
vian or non-Markovian depends largely on the nature of 
the initial state of the batbiiril. For baths at zero tem- 
perature, the decoherence is often non-Markovian. Ex- 
amples include the well known spin-boson models, and 
models of interacting spin environments that we already 
mentioned abov o^^i^^ . It has also been found in the case 
of a qubit interacting via the Fermi contact hyperfine 
coupling with a bath of polarized nuclear spinsi^. How- 



ever, for baths in thermal equilibrium i.e., at finite tem- 
peratures, the decoherence in the weak coupling limit is 
expected to be Markovian as in the case of the spin- 
boson model. This aspect was also seen in the case of 
interacting spin baths at finite temperatures^^"—. The 
associated Markovian decay rate was also found to in- 
crease with temperature and, unlike the case of bosons, 
was found to saturate at high enough temperatures. 

Since interactions between spins in condensed matter 
systems arc known to generate a whole range of com- 
plex thermodynamic behaviours, it is interesting to ask 
whether interactions between the bath spins generate 
novel behaviour for the coherence of a qubit coupled to 
such a bath. In this general context, various questions 
arise naturally: does the decoherence contain clear sig- 
nals about the underlying thermodynamic phase of the 
bath? Is Markovian decay to be expected for all interact- 
ing spin environments at finite temperatures, where the 
bath is in a statistical mixture at temperature T, and 
not in a pure stateiiii^ii^? And can decoherence or oth- 
ers physical quantities asociated with the qubit be used 
as a sensitive probe of the dynamics of the underlying 
spin environment? 

Here, we shall explore some of these issues by consid- 
ering the problem of the decoherence of a qubit induced 
by weak coupling to a model spin-bath with non-trivial 
dynamics arising from strong frustration and prepared in 
an initial state which is in thermal equilibrium. Highly 
frustrated or disordered baths were partially explored in 
lOllllI where complete predictions for the decoherence 



could be made only for the case of a one dimensional 
Ising bath or for spin shards in the infinite tempera- 
ture limit. Though these works illustrate the potential 
richness of interacting environments they also highlight 
the difficulty and limitations of using analytical methods 



to study these problems. We adopt a disordered spin 
bath described by the meanfield random Ising model in 
a transverse magnetic field. This environment is charac- 
terized by a rich thermodynamic phase diagram which 
shows a spin-glass to paramagnetic transition as a func- 
tion of temperature and external magnetic field^^. The 
total Hamiltonian for the system is given by 
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where Hb is the bath hamiltonian and Hsb denotes 
the spin-bath interaction. The Pauli matrices cTc and 
(Ji denote the spin of the qubit and the i-th spin of 
the bath respectively. The exchange energies J^ are 
quenched random variables with a probability distribu- 
tion Pi J) ~ ^=L— e^'-'"/^^', where the standard devia- 
tion of the distribution A sets the unit of energy. The 
qubit couples to the bath operator V" = ^^ ^^f, with 
a representing either the z or the x component of the 
spins. The coupling constants A are also taken to be 
quenched random variables. The physics of the bath 
is described by Hb, the random transverse Ising model 
(RTIM) which is known to have the following phase dia- 
gram: at zero temperature, for transverse fields smaller 
than a critical field h^ i.e., h < h^ the system is in a 
spin glass phased. For h > he the system is magne- 
tized with a gap in its spectrum. At finite temperatures, 
the spin glass order survives up to a critical temperature 
Tsg{h)^'^. A physical system which is expected to be 
reasonably described by the above model is LiHoxYi_xF4, 
where the Ho concentration, x, tunes the system between 
different physical regimes^!. It has been proposed as a 
minimal toy model to investigate the effects of quantum 
entanglement of spins^^. The critical temperature Tgg, 
that sets the value of the coefficient A, is of the order of 
lK2i, in the compounds. Consequently, the effective high 
temperature paramagnetic phases we study in this paper 
correspond to rather low real temperatures, where we 
expect the spin environment to dominate and phononic 
contributions to the decoherence can be neglected. There 
have also been propositions for engineering model mag- 
netic environments with very small energy scales using 
cold atomsH. 

We note that in the spin glass phase, only the correla- 
tions involving the a;-component of the spin exhibit spin 
glass features whereas the correlation of the Sz compo- 
nents merely exhibit a gap. Clearly, it would be inter- 
esting to explore the effect of these phases on the de- 
coherence of the qubit. Moreover, due to the inherent 



anisotropy of the bath, we expect the decoherence to be 
dependent on whether the qubit couples to the x or z 
components of the bath spins. In this paper, wc ad- 
dress these questions in the limit of weak qubit-spinbath 
coupling. Since this problem cannot be studied analyti- 
cally, we use numerical exact diagonalization methods to 
calculate the bath cigenstates and eigenvalues and con- 
sequently, the resulting decoherence. We obtain a rich 
spectrum of results for the coherence of the qubit and in 
particular, a non-Markovian i.e., a non-exponential decay 
of the coherence at finite temperatures. 

The paper is organized as follows: in SecHH we present 
the weak coupling formalism used to calculate the deco- 
herence of the qubit. We present our numerical method 
in Sec. IIIII and use it to study the analytically tractable 
case of decoherence induced by an Ising chain spin-bath 
so as to benchmark our method. We compare our results 
for the Markovian decoherence rate with known analyti- 
cal results for an infinitely long chain, so to establish the 
importance of the finite size effects introduced by our nu- 
merical method. In Sec. IIV Al we present our results for 
the decoherence in the weak coupling regime for the long 
range Sherrington-Kirkpatrick model which has an Ising 
spin glass phase at low enough temperatures and present 
a comparison of our results with some analytical results 
which are known in the high temperature paramagnetic 
phase. In Sec. IIV Bl we present in detail our main results 
for the decoherence induced by the RTIM spin-bath for 
two different couplings of the qubit to the bath. We con- 
clude with a discussion of our results in Sec. IVl 



II. WEAK COUPLING FORMALISM 



In this section, we summarize the weak coupling 
formalisntii used to calculate the decoherence. This ap- 
proach is valid provided the energy associated with the 
qubit-bath coupling is smaller than all the scales of the 
bath. We first assume that at time t = 0, the combined 
system of the central spin and the bath have a factoriz- 
able initial density matrix: fl = p{0) ® ps- Depending 
on the coupling of the central spin to the bath operator 
y , the qubit is in a pure state whose basis vectors are 
defined by 1^}") ^ a\ ^'^) + /3| ^'') (p(0) = |i/^")(V'1) 
where I ^^) = I ^), I ^'=> = I ^> and | ^^) = | t), 
I ^■^) = I i). The bath is chosen to be at thermal equi- 
librium with temperature T = 1//3 leading to a density 
matrix 



PB = 



-PHb 

z 



(2) 



where Z — Tv cxp(—f3HB) is the bath partition function. 
The time evolved reduced density matrix is given by 



pit) = lap I ^-){^- I + l/^n ^-) (^- I + M(i)a*/3| ^'^) (^"^ | + A/(t)*a/3* | ^'^) (^"^ | 



(3) 



r 



where the factor 



A'f(i)=Trfe-*(^^+^°)*PBe' 



j(HB-y")t 



(4) 



is a measure of the decoherence induced by the bath at 
time A^iii. Note that Tr denotes the usual trace as Hb 
and y" arc operators in the bath Hilbert space. For 
weak coupUng to the environment i.e., for smaU A, we can 
use the super-operator formahsmi to obtain the following 
form for the decoherenceii, being the Laplace transform 
of M{t) 



M{z) 



-I I dte"^M(t) 
'o 



(5) 



where ^ is a complex variable with Imz > 0. As shown 
in |ll| . this Laplace transform can be written as 

M(z) = [z-S](z)]-^ (6) 

where the self-energy S up to second order is given by 
E2(^) = 2Tr(y>B)- 

2i / dte"* [(l^"(t)y")e + (VVit))^] (7) 

"'0 

where the connected correlation functions are defined as 



(i/"(t)y")c = {V'Vit)) - {v''{t)){v'') 



(8) 



The coherence M can thus be written in terms of the real 
functions A2 and r2 defined by 



A2{E)-iT2iE) 



lim S]2(£: + iry) (9) 



where E is real. This then leads to the weak coupling 
result 



'its 



e(t)M(t) = — / dE -^- -- 

"■ ^ "■ ^ 271 J E-h2{E)+iT2{E) 



(10) 



where 0(i) is the Heaviside step function and A2 and r2 
satisfy standard Kramers-Kronig relations. In general, 
r2(-E) might have analytic as well as non-analytic parts. 
It is interesting to note that in the weak coupling limit, 
the decoherence is essentially dictated by r2 [E) which as 
can be seen from ([7]) is proportional to the symmetrized 
dynamic structure factor of the bath. 

Typically for spin baths, we encounter situations 
where, one can have a net magnetization along a cer- 
tain spin direction {V"") = k 7^ 0, and/or magnetic or- 
dering where {V°'{t)V°')c = M + /(O with n and /i being 



constants independent of time, and f{t) a time-varying 
function. The resulting self energy in the Laplacian vari- 
ables then takes the form 



M 



£2(2) ^2K + A-+g{z) 



(11) 



where g{z) is some analytic function of z. For the asymp- 
totic decoherence, the first two terms generate oscilla- 
tions resulting in a decoherence of the form 

M{t) = cxp {2iKt) cos {2y/jlt)M{t) (12) 



wher( 



and 



.11 



lnM(.).-^/.i.^^i(f^r2(i.) 



T2{E) = [f{E) + f{-E)] 



(13) 



(14) 



being / the Fourier transform of the function f{t). ((T3)) is 
also known as the time-convolutionless projection opera- 
tor (TCL) approximatioii^^. The true intermediate time 
decoherence is expected to be lightly modified with re- 
spect to the result predicted by p^ . but it is expected to 
be qualitatively similar to that predicted by (fTO|) . which 
is valid for all t. In the absence of any intrinsic dy- 
namics of the central spin, the second order aproxima- 
tion used here, leads to an equation for the decoher- 
ence which resembles that obtained for the spin-boson 
model^l. Though PH)) is exact for the spin-boson case, 
here it is valid only in the limit of weak qubit-spinbath 
coupling. Morevoer, not all spectral densities r2(£') ob- 
tained from interacting systems can be simulated by the 
usual non-interacting boson baths. 

For the disordered systems studied here, the correla- 
tion functions (|S]) and ([7]) need to be averaged over the 
probability distributions of both the exchange interac- 
tions Jij and the coupling constants A,;. In the rest 
of the paper, the coupling constants A^ in ([T]) arc cho- 
sen to have the following disorder averages: A^ = and 
XiXj = {}? /N)6ij where TV is the total number of spins 
in the bath. This choice leads to the vanishing of the 
first order correction to the self energy since the disorder 
average of k is zero. The disorder averaged r2, is then 
directly r elated to the connected local spin correlation 
functions {V°-V°-)c, and it is evaluated using the spectral 
representation for the dynamical structure factor 



gaa^^^ 



dtcxp{iEt){V°-{t)V°-)c 



M N 2" 



27rA2 

NM ^ Z( 

771—1 

2TTfiS{E) + f{E) 



(15) 



where M is the number of reaUzations of disorder for 
the exchange interactions Jij and tti^ is the disordered 
average of the square of the local magnetic moment. Our 
method for evaluating ([T5|) is presented below. 



III. NUMERICAL METHOD 

To compute T2 for the RTIM spin-bath we use ex- 
act diagonalization methods. The advantage of the full 
diagonalization is that it permits us to calculate finite 
temperature quantities as well. We first study the coher- 
ence of the qubit coupled to a random Ising chain bath 
described Hb with the external field set to /i = , where 
known analytical results can be used to benchmark the 
implemented numerical method. The interactions in the 
bath are confined to nearest neighbors on the chain and 
the spin-bath coupling is taken to be in the a = z direc- 
tion. 

We consider baths with the number of spins A^ varying 
from 3 to 10 spins. For a given set of exchange interac- 
tions, the resulting hamiltonian which is a 2^ x 2^ matrix 
is exactly diagonalized and the eigenenergies and func- 
tions are obtained. These are then used to compute the 
dynamical structure factor and T2{E). Since the interac- 
tions are quenched random variables, we repeat the above 
procedure for several thousands of realization of disorder 
and we obtain the disorder averaged T2. In Fig. [U we 
plot r2{E) for the case of the Ising bath for different 
sizes of the bath, and fixed temperature T = 0.01 A. The 
analytical results for the thermodynamic case are also 
shown in the graph. 

Note that as the number of spins N in the chain in- 
creases, the numerical results converge to the analytical 
curve obtained for a thermodynamic bath. At high E, 
r2 {E) is almost independent of the size of the bath. This 
is due to the fact that the typical energy bonds are of or- 
der Jij {0{A)) (independent of the size of the bath), 
which arc sufficiently small to produce any modification 
in the high energy spectrum. At low energies, size effects 
are most significant, since at finite sizes there are no ex- 
act cancellations that lead to the vanishing of r2(£') at 
E = 0. Nevertheless, the spectra shows good uniform 
convergence as N increases. At higher T, when the ther- 
mal energy is higher than the typical energy of the bonds, 
the convergence is even better and T2{E) is roughly in- 
dependent of the size of the bath. In Fig. [21 we plot the 




FIG. 1: (Color online) A\-'^r2{E) as a function of E for 
T = O.OIA and different A''. Note that as N increases, the 
numerical results approach the analytical result (An.). A sim- 
ilar result is found for the SK long range system. 




FIC. 2: AA"^r2(£) as a function of E for different T for 



variation of F2(£') for different temperatures. As shown 
in Ref. Illl it is the low frequency part that is mainly 
affected by the temperature effects. 

At high T, r2{E) is basically a broad peak of width 
A, centered a.t E = 0. This is due to the thermal fluc- 
tuations which result in uncorrelated spins, and T2{E) 
merely reflects the distribution of the bonds Jij and the 
structure of the energy levels. Both temperature and size 
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FIG. 3: AA"^r2(0) as a function of 1/N'^. Note the excellent 
agreement of the extrapolated value of r2(0)and the analyti- 
cal one, except at low T {T < 0.2A). 
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FIG. 4: AA^'^r2(0) (extrapolated to the thermodynamic 
limit) as a function of T. Good agreement with analytical 
results (dotted curve) except at T < 0.2A (inset). 



effects are seen to have an impact on the low frequency 
part of r2 {E) with the consequence that the asymptotic 
decoherence is more sensitive to thermal fluctuations and 
finite size effects. 

We use our numerical results for T2{E) to calculate 
the decoherence M{t) given by ([13]). A comparison of our 
results with the analytical results obtained for the infinite 
chaii*ii permits us to gauge the importance of finite size 
effects. As for the infinite chain case of Ref. [ill, the 
decoherence for the finite chain can also be analyzed in 
terms of three regimes described by the two characteristic 
times: t^ = A"-'^ and tp = fi = T~^ . At short times, the 
coherence is given by the universal gaussian. In M(t) = 
-<2(27r)-i / dET2{E), given by the sum rule of r2iE) 



dE r2{E) =47rAVA 



(16) 



The asymptotic Markovian regime, where the coher- 
ence decays exponentially is characterized by the values 
of T2{E) at E ^ (for t -^ oo, \nM{t) ex -r2(0)t). 
At high T, M{t) remains practically constant up to 
t ~iA, where a final Markovian regime arises. When 
T < A, the gaussian regime is followed by an intermedi- 
ate time regime where the coherence decay as a power 
law, before reaching the Markovian regime. A care- 
ful inspection of our numerical data shows that, r2(0) 
scales as 1/iV^ (Fig. [3]). The extrapolated values of 
r2(0) for A^ — > OD as a function of T arc plotted in 
Fig. [31 where it can be seen that there is an excel- 
lent agreement between the analytical values given hy^ 
r2(0) = 27rA2/dJF(J)2(l-tanh(/3J)2), and the ex- 
trapolated values except for T ^ 0.2 A (inset of Fig. [4]). 
At very small T, since both size effects and temperature 
effects are non negligible at low energies, we surmise that 
the accessible system sizes are not sufficient to reach the 
scaling regime. 

To summarize, we see that our numerical method does 
reproduce the expected analytical behavior for the spin 



chain bath except at low temperatures where, the Marko- 
vian rate is overestimated in our approach. Nonethe- 
less, our results reproduce the existence of the different 
regimes. We would like to draw attention to the fact that 
the benchmark case of the chain is in a sense a worst case 
scenario. In fact, for a 1-d spin-bath one expects to have 
stronger finite size effects compared to the models to be 
studied in the rest of the paper: i.e., infinite-range mod- 
els which are known to have much reduced finite size 
effect922^. 

Before leaving this numerical method section we would 
like to note that our method differs substantially from 
others that are based on the solution of finite size clus- 
ters and cluster expansions. Those methods usually deal 
with a non-disordered model hamiltonian, which is solved 
for a small system size. The finite size effects are dealt 
with by implementing various cluster expansion schemes. 
In those cases, the finite nature of the cluster provides a 
discrete pole structure with a finite low energy cut-off 
and hence a finite recursion time that provides a long 
time cut-off for the study of the decoherence effects. In 
our case the situation is different. The effect of the bath 
enters through the spin susceptibility, and this quantity 
in a disordered model is computed through the disorder 
average. Since we have a continuous (Gaussian) distri- 
bution of couplings, the disorder-averaged susceptibility, 
though computed on finite size clusters, does not have 
a finite low-energy cut-off. Therefore, the decoherence 
does not have a finite recursion time. All the systematic 
errors in our approach are due to the second order ap- 
proximation, which is safe so long the coupling between 
the qubit and the bath is small, and to the finite size of 
the clusters that are diagonalized to compute the suscep- 
tibility of the bath's Hamiltonian. This latter quantity 
can be, nevertheless, extrapolated to large system sizes, 
as we shall see latter. This is validated by the benchmark 
of the method against an analytically solvable case that 
we describe in the beginning of next section. 



IV. MODELS OF SPINS BATHS 

In this section, wc apply our numerical method to 
study the decoherence induced by interacting quantum 
spin-baths which do not have conventional magnetic or- 
der like ferromagnctism or antiferromagnetism. The 
models studied all have infinite range interactions: i.e., 
all spins interact with each other which leads to the pres- 
ence of a significant geometric frustration. Similar mod- 
els have been studied in Rcf. \1Q, where the focus was on 
the nature of the distribution of the energy level spacings 
and their impact on the decoherence. A standard lore is 
that thermal fluctuations in the weak coupling regime 
always leads to a Markovian i.e., exponential decoher- 
ence. Here, we will discuss cases, where the decoherence 
remains highly non-Markovian in certain finite tempera- 
ture paramagnets. 
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FIG. 5: (Color online) AA"^r2(£;) for different values of T 
and N = 9, for the SK bath. The results are qualitatively 
similar to that of the chain. 



A. Infinite-ranged Ising batii witii h — 

We study the decoherence induced by the infinite- 
range random Ising bath also known as the Sherrington- 
Kirkpatrick (SK) model^^. For a gaussian distribution 
with zero mean of the exchange interactions between the 
spins, the system is paramagnetic except below the spin- 
glass transition temperature T < T^g = A where the 
system develops spin glass order. Note that for the in- 
finite range model, extensivencss of the free energy re- 
quires as to scale the interactions J^- —J- Jij/N . In this 
case, there is no first order contribution to the self en- 
ergy and the second order contribution is given directly 
by the full correlation function as opposed to the con- 
nected correlation function of ([5]) i.e., both k = ^ = 0. 
Moreover, r2 is proportional to the probability distribu- 
tion of the local magnetic fields as discussed in Ref. |Tl|. 
This distribution is straightforward to evaluate analyti- 
cally in the paramagnetic phase, but numerical methods 
are required to obtain the same in the spin glass phase. 
Earlier work only discussed the evolution of the Marko- 
vian rate close to the spin glass transition temperatursii. 
Here, we compute the coherence for all times and find an 
asymptotic Markovian regime at all finite temperatures. 
Our results are shown in Fig. [S] and we see that they 
qualitatively resemble the results obtained for the Ising 
chain. Moreover, our method reproduces known analyti- 
cal results in the paramagnetic phase and also approaches 
the expected linear behavior of F at low frequencies as 

The numerical results for In M{t) are shown in figurc[Sl 
At short times, t ^ t/\ = 1/A, the sum rule obeyed by F2 
ensures that M{t) decays as a gaussian akin to the free 
spin bath. For times t > t/^, temperature starts playing 
a relevant role. At high T > Tgg, since ^2{E) is a peak of 
width 0(A) cf. Fig. [S] M{t) decays as a Markovian for 
all times t ^ tA- For T <^ Tsg, the asymptotic Marko- 
vian regime is preceded by a power law regime resulting 
from the partial linear behavior of T2{E) for small E. As 
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FIG. 6: (Color online) (A/A)Mn M(i) for the SK bath and 
the comparison with the free spin-bath (FSB); (a) for short 
times, where it can be seen that at T = A, there's a change 
of the curvature in In M(i); (b) for long times. 



T — > 0, this intermediate power law regime is expected 
to extend to the asymptotic regime. This is however, 
hard to infer from the numerics, since finite size effects 
smear the value of r2(0) and hence, the linear behavior 
of r2 (E) predicted by the TAP-method calculation [30[ . 
Nevertheless, we see that as T decreases from the para- 
magnetic phase, r2(0) decreases and the linear contri- 
bution to r2(£') become more significant. Extrapolating 
the results of r2(0) to the limit of A^ — >■ 00 to obtain 
the Markovian decay rate in the thermodynamic limit, 
we find that r2(0) scales 1/A^ cf. Fig. [7] (instead of the 
1/A^^ scaling in the Ising chain). The resulting Marko- 
vian rate is plotted as a function of T in Fig. |S1 and 
perfectly matches the analytical values predicted in [11| 
in the PM-phase. As before, there is a mismatch of the 
results in the very low temperature regime, probably be- 
cause our system sizes do not access the scaling regime. 




FIG. 7: Scaling plot for AA~2r2(0) as a function of l/N for 
the SK bath. 
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FIG. 8: AA~ r2(0) (extrapolated to the thermodynamic 
limit) as a function of T for the SK bath. The dashed curve 
represents the analytical results in the paramagnetic phase. 



As in the case of the Ising chain, we find that interac- 
tions between the spins do lead to longer coherence times 
as opposed to the case of a free spin bath. The only visi- 
ble effect of the spin glass phase transition is the point of 
inflexion in the curvature of In M(f) for f ~ tA when the 
asymptotic Markovian regime takes over. This curvature 
is negative for T > T^g and positive for T < T^g. The 
absence of a radical change in the decohcrcncc is not sur- 
prising since the qubit does not couple to the operator 
which corresponds to the spin glass order parameter. 

In the following section, we study the RTIM which 
is known to have both, a quantum phase transition at 
T = 0, and a finite temperature classical phase transi- 
tion. We study two different couplings between the qubit 
and the bath spin operators and show that the nature of 
these couplings has highly non-trivial consequences for 
the dccoherence. 



B. Sherrington-Kirkpatrick model in a transverse 
field ft / 

The Sherrington-Kirkpatrick model in a transverse 
field (RTIM) is a more interesting case than the previ- 
ously studied ft, = case. At T = 0, the system un- 
dergoes a phase transition from a spin glass state for 
h < he — 1.44A to a gapped phase for h > he- As 
temperature is increased, the spin glass phase disappears 
at a finite temperature which depends on the value of 
the magnetic field. The magnetic field ft is therefore, ex- 
pected to deeply influence the behavior of r2{E), and 
hence the coherence of the central spin. Since the model 
is not isotropic, it is important to note that the spin glass 
order exists only along the x component of the spin. To 
probe the physical ramifications of this anisotropy, we 
examine two different couplings of the central spin to the 
bath: a coupling of the spin operators in the direction 
of the field (a ~ z) and a coupling perpendicular to the 
field [a = x). As we will show below, these lead to radi- 
cally different predictions for the decoherence. Previous 
numerical studies^ of this model have shown finite size 
effects to be rather minimal. Given the numerical com- 
plexities in the vicinity of the phase transition, we limit 
the scope of the present study to consider two values of 
the field ft = O.lftc and ft ~ 2ftc, that set the system in 
qualitatively different regimes. 



1. Coupling parallel to the field (case a = z) 

Since a magnetic moment is present in the z direction 
for ft 7^ 0, the T2{E) is given by the full connected cor- 
relation function ([5]). As shown in (|15p . this contributes 
a singular term ex IJ,S{E) to ^2{E) which then leads to 
oscillations in M{t) (equation [T^]). An unambiguous way 
to extract the coefficient /i is via the sum rule satisfied 
by the dynamical structure factor—. This, however, is 
numerically cumbersome and in the rest of the paper, 
we do not deal explicitly with these singular terms, since 
they only induce oscillations and instead concentrate on 
the non singular part of ^2{E) that leads to a decay of 
the coherence. The T = behavior of r2(£') when a 
small field ft= O.lftc is applied is shown in figure IHl We 
note that even for such small fields, T2{E) is radically 
different from the earlier case of ft = 0. 

Our results for r2 (E) obtained for A^ = 9 are shown in 
Figs. [^ and [TUl for the small and large field case respec- 
tively. In the former, with ft = O.lftc, r2(£') exhibits a 
gap of order 2ft at T = (numerically, the T = curve 
is indistinguishable from the T ~ 0.1 A curve shown in 
Fig. [S]). At high T, T2{E) shows a pronounced peak 
around E ~ ±0.1ft and broader peaks of width sa 4A in 
the background. Similar features are seen at high mag- 
netic fields (Fig. [To)) : existence of a gap 2ft at T = 
and a two peak structure of width 4A at high tempera- 
tures. For high enough T, the susceptibility at low fre- 
quencies becomes linear with a temperature dependent 




FIG. 9: (Color online) AA"^r2(£) for diflterent T, with N = 
9 and magnetic field set to h= O.lhc- Note the completely 
different behavior from the case h = (Fig- O at small T, 
r2(-E) is gapped (right inset, r2(-E) in a logarithmical scale). 
For higher T (left inset), r2(-E) is linear at low energies. We 
note that the T = curve is indistinguishable from the one 
at r = 0.1A. 




FIG. 10: (Color onhne) AA~^r2(£;) for different T, with fixed 
A*' = 9 and field h= 2/ic(the T = curve is not plotted here 
since is indistinguishable from that of T = O.IA). Right inset: 
behavior for low E. Left inset: high frequency peaks. The 
missing spectral weight is carried by the 5{E) term described 
in the text. 



slope, r{E) = a{T)\E\ as E ^ 0. Note however, that in 
the high field case the S{E) tcrni arising from moment 
formation carries most of the spectral weight at T — 0, 
as the spectral weight of the regular part of T2{E) de- 
creases dramatically. The principal difference between 
the high and low field case is the presence of the sharp 
peak around E ^ 2h in the low field case. We have stud- 
ied the finite size scaling of T2{E) and we find that the 
features mentioned above are robust to finite size effects. 
Moreover, contrary to the previous cases, here we obtain 
a good convergence with system size in the low frequency 
regime. 

Interestingly, our results bely the naive expectation 
that the behavior in the high temperature paramagnetic 
phase be independent of the values of h considered here. 
This can be attributed to the fundamental difference be- 
tween the structure of eigenfunctions and eigenvalues in 
the two cases: for h = O.lhc, they correspond to that of a 
spingiass system, while the high field case essentially cor- 
respond to spins in a strong field where the interactions 
between the spins can be viewed as a perturbation. At 
high enough temperatures, all eigenfunctions and eigen- 
values contribute equally to the r2(£') thus the large T 
regime actually reveals the "geometric" underlying struc- 
ture of the Hamiltonian, rather than a naively expected 
universal paramagnetic regime of free spins. This obser- 
vation remains relevant for all the different spin-models 
analyzed in this work. 

The results for r2(£') indicate a very rich evolution 
of the coherence of the central spin, particularly at long 
times. We first analyze the results for M{t) plotted in 
Figs. [H (high T) and [12] (low T), for h = O.lhc- As 
before, for t ^t^, the sum rule obeyed by T2{E) P^ 
dictates an universal gaussian decay for the decoherence, 
independently of h, A and T. This is indeed verified by 
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FIG. 11: (Color online) {A/\flnM{t) with h = 0.1/ic. The 
upper graph (linear scale) , shows an intermediate time Marko- 
vian regime { t > t ':^ tm), while the lower graph (logarithmic 
scale) shows a final power law regime t > t c^i tp). Insets: 
In M{t) for lowT. 



our results. 

For r = (insets of Fig. \TT\ , the decoherence is only 
partial due to the presence of a gap in r2 (E) and shows 
oscillations with a frequency proportional to the size of 
the gap. Similar oscillations are also seen at very low 
temperatures. The short time regime t < t/^ is followed 
by a power-law regime for t > tp which saturates to a 
finite value at long times. As T increases, the short time 
and asymptotic power law regime are separated by an 
intermediate Markovian regime for i™ < i <C ip. For 
comparison, the coherence for a bath with A^ = 10 spins 
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FIG. 12: (Color online) {A/\flnM{t) with h = O.lhc The 
upper graph (linear scale) , in contrary to the high T coherence 
of Fig. 111! does not show an intermediate Markovian regime, 
while the lower graph (logarithmic scale) shows a final power 
law regime in which tp increases with T. 
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FIG. 13: (Color online) Coherence of the central spin as a 
function of time, for a bath of A^ = 10 spins for h = 0.1/ic and 
a bath-qubit coupling A — 0.2A. The free spin case (FSB) is 
also plotted for reference. 



is plotted in figure [T3] various T. As in the case for h = 0, 
not dramatic change is seen as one traverses the spin glass 
transition temperature. The existence of an asymptotic 
power law regime at finite temperatures completely defies 
the conventional lore that the asymptotic decoherence 
at finite temperatures is Markoviani^iiS. A systematic 
analysis of our numerical results indicates that for low 



fields, tp is an increasing function of temperature. This 
implies that in the high temperature paramagnetic phase, 
the power law regime is pushed to ultra long times and 
the decoherence is Markovian for realistic times. Though 
the exact values of tp and t,„ are subject to finite size 
effects, we have studied various system sizes and find 
that the intermediate Markovian and power law regimes 
will survive in the thermodynamic limit. 

Before moving on to the high magnetic field case, we 
shall present an intuitive physical picture of the very low 



field case. We first recall the physical picture of the spin- 
glass ground state at T = and h = 0, that was discussed 
in the numerical investigation of Ref. [22. There it was 
argued that the spin system can be viewed as a coexisting 
collection of large non-frustrated clusters of spins with 
few remnant strongly frustrated "dangling" spins. Thus, 
the dangling spins experience a distribution of effective 
magnetic fields h^ff which, from the work of Sherrington 
and Kirkpatrick is known to have a linear distribution. 
This linear is directly reflected in the r2{E) ^ \E\ that 
we discussed before. 

Now, when a small external magnetic field is turned 
on, the unfrustrated spins of the clusters will remain es- 
sentially unaffected. In contrast, the dangling spins will 
align in the direction of max(ft,off, /i). This means that 
the finite external field h will act as a low frequency cut- 
off, pushing spectral weight towards higher frequencies 
and thus opening an /i-controUcd gap in T2{E) around 
E = (see the low-T curve of the right inset of Fig. [9]) . 
This behavior is in stark contrast with the reference case 
h = presented in section IIV Al where T2 (E) remains 
always gapless (Fig. [5]). The dramatic change in the low 
E regime implies also the qualitative change in the long 
time behavior of the decoherence that we described in 
the previous section. 

This line of argumentation also allows us to qualita- 
tively understand the origin of the linearity of T2 {E) at 
high T. In this limit, from ([T5|) we see that all eigenval- 
ues and eigenfunctions contribute to r2 {E) . On the other 
hand, when /i = 0, the data of Fig. [5] shows that T2{E) 
is finite as £' — > 0, so there are plenty of contributions 
at arbitrary low frequencies, which originate in pairs of 
quasi-degenerate eigenstates \n) and \m) with eigenener- 
gies En ~ Em, which also have a non- vanishing matrix 
element s^m = {"'T-Wil'i^^)- When a small external field ( h 
<C /ic ) is turned on we may consider it as a perturbation 
an perform the following qualitative analysis: The (2x2) 
block of H in the subspace of n and m will read: 



^nni 



Eo 



Eo 



where Eo is the quasi-degenerate energy of the states. 
The effect of the small h is to lift the degeneracy of the 
pair of levels resulting in a transfer of spectral weight 
from E Ri to higher frequencies. Extending the analy- 
sis to all pairs of quasi-degenerate states contributing to 
T2{E) implies the immediate collapse of the finite value 
of r2{E = 0) to zero when hj^ 0, as is seen in our results 
of Fig. [TH Moreover, the new eigenvalues of Hnm are 
sa ±hs^j^^, thus with the reasonable assumption of a fea- 
tureless distribution of the s^^ matrix elements, we con- 
clude that the degeneracy lifting is uniformly distributed 
and directly proportional to h. This implies that the be- 
havior of r2{E) at very low frequencies should be linear 
in E with a slope roughly given by 1/h. This analysis is 
in qualitative agreement with the low field data shown in 

Fig. [13 

At high magnetic fields (Fig. [T5]) . the coherence 
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FIG. 14: (Color online) AA"^r2(£) as a function of E for 
different values of the field h with TV = 9 and T = lOA. 
Inset; linear behavior of r2(£') for low energies. 
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FIG. 16: The exponent e {M{t) « i '^) extrapolated to the 
thermodynamic limit, as a function of T for h = 2hc. 



2. Coupling transverse to the field (case a — x) 
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FIG. 15: (Color online) {A/ \f In M{t) as a function of time 
for h= 2hc. For low T, since r2{E) is gaped, the decoherence 
is very weak and shows oscillations (inset). As T increases, 
the coherence decays asymptotically as a power law. 



presents a different qualitative behavior. At low T, since 
T2{E) is nearly gapped, and most of the spectral weight 
is in the S{E) term, the central spin decoherence is very 
weak and shows oscillations. At high T, in contrast to 
the low h behavior, there is no well defined intermedi- 
ate Markovian regime, but we recover a power-law decay 
{M{t) DC t~^) at asymptotic times. This coefficient e, 
which characterizes the power law decay, increases with 
T and saturates to a finite value at very high tempera- 
tures, proving that the qubit decoheres faster as temper- 
ature increases. The value of the exponent e extrapolated 
to the thermodynamic limit, as a function of T, is shown 
in figure 1161 



We now explore the case where the qubit couples di- 
rectly to the bath operator which is related to the spin 
glass order parameter. In the spin-glass part of the phase 

diagram {h < he), we now have {a^) = and (cr^) ^ 
Vi. As in Sec. IIVB 11 T2{E) is again a sum of a regu- 
lar and a singular contribution. The singular part is a 
S{E) contribution whose strength is related to the spin 
glass order parameter or Edwards Anderson parameter 
<ZEA, rNR = (47rgEA-E^MF)^(^)- The non-regular 
part induces oscillations in M{t) and the regular part re- 
sults in the decay of M(i) ([T^. The estimation of ^ea 
is a central problem in spin-glasses and as expected our 
numerics do not produce well converged results of the 
values of this parameter, qea is typically estimated by 
systematically studying T2{E) for different sizes of the 
system, as in Ref. [22- Though it is reasonably straight- 
forward to estimate qea at T = 0, this is not the case at 
finite temperatures where the delta peak is broadened by 
both finite size effects and thermal excitationsi^ 

When h — 0, the operator V^ is constant of motion 
since [Hb^V^] = resulting in a temperature indepen- 
dent r(^) = AnX^A-^SiE). As discussed in Sec. IH 
this results in a purely oscillatory behavior of the form 
M{t) = cos(2AA~^i).When the field h is turned on, spin 
flips arc allowed, [Hb,V^] ^ and T2{E) has a richer 
structure. For small h (h <ti he) and T = 0, in addition 
to a S{E) term r2{E) has a regular part r^(^) oc \E\ 
as £' — s- and has a peak of width 4A, centered around 
E = 2hS^. As temperature is increased, since qea is ex- 
pected to decrease, a part of the spectral weight of the 
6{E) term is transferred to the regular part F^ which is 
now nonlinear at low energies (Fig. I17(). 

We can now construct a picture of the coherence in the 
spin- glass phase based on the above description of F2(£'). 
The short time behavior is governed by a gaussian. Due 
to the afore-mentioned problem with the broadening of 
the delta term, we have not been able to obtain clear pre- 
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FIG. 17: (Color online) AA-^r2(£') for N=9 and field h= 
0.1/icfor different T. For T — 0, r2(-E) increase linearly for 
low E and has a maximum at E = 2h. As T increases, spectral 
weight from the 5{E) (discussed in text) is transferred to low 
energies. Inset: V2[E) for h — 2hc showing thermal smearing 
of the gap. 
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FIG. 18: (Color onhne) AX-'^\n M{t) for N=9 and field h= 
2hc for different T. For T = 0, since r2(ii') is gapped, the 
decoherence is partial but asymptotic Markovian decoherence 
is recovered at finite T. 



dictions for the intermediate and asymptotic decoherence 
for low h and low T. This is intimately linked to the un- 
derlying spin-glass order since in this case some of the 
contribution to the singular part exists only in the ther- 
modynamic limit, as opposed to the preceding section, 
where the singular term arises from a straightforward 
magnetization of the underlying system and hence one 
did not have to deal with broadening induced by finite 
size. Nonetheless, for T = and /i <C /ic it is possible 



to have a qualitative picture of the coherence. As in the 
zero-field case, the linearity of T2 (E) a.s E ^ leads to 
a power-law asymptotic decay of the coherence. At fi- 
nite temperatures, we expect the asymptotic decay to be 
Markovian. 

For high h {h > he), there is no spin glass order and 
hence no singular contribution to T2{E). At T = 0, the 
function T2{E) is gaped, and it is centered in E — ±2hc, 
with a width of ~ 4A. As T is increased, thermal ex- 
citations appear within the gap (Fig. [T7)) . Thus, in the 
former case the decoherence is partial (Fig. [T8)) and we 
see oscillations with a frequency given by the gap size. 
As T increases, thermal excitations start to fill the gap 
and M(t) decays in a Markovian way given by the value 
of r2(0). This behavior conforms to the usual expecta- 
tions of asymptotic Markovian decay at finite tempera- 
tures and is very different from the small h case. 



V. CONCLUSIONS 

In this paper, we have used exact diagonalization 
methods to study the decoherence of a central spin in- 
duced by spinbaths with random interactions. We find 
that the asymptotic decoherence is intricately linked with 
the nature of the interactions in the bath. Moreover, for 
a given set of bath parameters, the decoherence crucially 
depends on the nature of the coupling between the qubit 
and the bath spins. For the cases of the random trans- 
verse Ising model studied here, we find that the under- 
lying nature of the eigenstates of the bath hamiltonian 
play a preponderant role in determining the decoherence 
even in the finite temperature paramagnetic phase. More 
precisely, we find that the decoherence in some of the 
finite temperature paramagnetic phases is strongly non- 
Markovian. Wc emphasize that standard Markovian ap- 
proximations used to obtain the density matrix and the 
decoherence at finite temperatures should not be used 
blindly as they can lead to highly misleading results as in 
the present case of disordered interacting systems. Unfor- 
tunately, we have not been able to study the decoherence 
in the vicinity of the phase transitions in the bath ( both 
at zero and finite temperatures) . This requires the use of 
other numerical and analytical methods which is beyond 
the scope of the present work. It will be interesting to 
study if these features are seen in other highly frustrated 
spin baths and whether this non-Markovian behavior sur- 
vives when order corrections to the self-energy are taken 
into account. These questions are left for future work. 
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